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1 Introduction 


CP violation in neutral meson-antimeson systems arises from three contributions: 
direct CPV, which is CP violation in a specific decay, and indirect CPV from mixing 
and from interference of decays with and without mixing. 

Mixing in decays of to CP eigenstates {K~^K~, tt+tt”) results in a different effective 
lifetime than that of decays to flavour eigenstates (iC“7r+). The observable ycp is 
defined as 


_ r(L)° ^ K-71+) 
~ t{D^ K+K-) 


( 1 ) 


CP violation gives rise to a difference of effective lifetimes of and decays to 
the same CP eigenstate. The asymmetry Ar can be defined as 


r(L)° ^ K-K+) - t{D° K+K-) 

Ar = — =0 - • ( 2 ) 

t{D K-K+) + r(L)0 ^ K+K-) 


It follows that if no indirect CP violation is present, Ay becomes zero. 

The Standard Model prediction for CP violation in charm is of the order 10“^ [T]. 
Any measurements of a larger value would indicate signs of New Physics, making this 
field an interesting area of study. 


2 Time-dependent analysis of —)► K^K , tt+tt 

Through a measurement of ycp and Ay we measure indirect CP violation, which is 
common to all decay modes. The quantities ycp and Ay are extracted via proper 
decay time measurement. We present here the final result of the update of the same 
analysis, performed by the Belle Collaboration in 2007 [2]. The new results cover the 
entire Belle data sample, corresponding to 976 fb“^ collected at or near T(4S), T(1S), 
T(2S), T(3S) and T(5S) resonances. An improved analysis method is implemented, 
which recognises two different conhgurations of the Silicon Vertex Detector (SVDl 
and SVD2) that were used during the data taking |3], and accounts for the polar 
angle dependence. 

2.1 Event Selection 

The is required to come from a decay chain D*^ —The charge of the slow 
pion provides the necessary tag on the flavour of the charm meson. Additionally, it 
enables us to set a constraint on the total energy released in the decay, q = — 

m{D^) — which provides excellent background suppression. 

To select kaons and pions, standard Belle particle identification is applied [1]. The 
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daughters of the are rehtted to a common vertex. For the vertex £t of the D*^, 
and 7r+ are also fitted to the interaction point. In both cases, a confidence level 
greater that 10“^ is required. To exclude D^s coming from decays of B mesons, a 
cut on the CMS momentum of the D*~^ is imposed. It is required that Pcms{D*~^) > 
2.5 GeV (3.1 GeV for T(5S)). Gandidates for D^s are selected based on the mass 
m{D^) and q, with a window imposed around the nominal value. For m{D^), the 
allowed deviation is ±2.25 um, where um is the R.M.S. width of the peak (um) 
and depends on the decay mode and SVD conhguration. Typical values are aM ~ 6—8 
MeV. The window in q is ±0.66 MeV (0.82 MeV) for SVDl (SVD2). The selection 
criteria are optimised in order to achieve the minimal statistical error on ycp- The 
proper decay time of the is calculated from the projection of the vector joining the 
two vertices L onto the momentum vector: t = mnom{D^)L-p/p^, where mnom{D^) 
is the nominal mass. The proper decay time uncertainty is evaluated from the 
error matrices of the production and decay vertices. It is required at < 440 fs (370 
fs) for SVDl (SVD2) to reject candidates with poorly determined proper decay time. 
The proper decay time distribution is parameterised as 

F{t) = ^ J - t')dt' ± B{t) , (3) 

where r is the effective lifetime, N is the signal yield, R{t) is the resolution func¬ 
tion and B{t) is the background distribution. Background is hxed from a £t to the 
sideband distribution. The position of the sidebands is optimised as to minimise the 
systematic uncertainty. The resolution function is constructed using a normalised 
distribution of at, combining hts for different bins. An additional offset parameter is 
introduced to correct for a certain misalignment of the SVD detector. As this param¬ 
eter is a function of the cosine of the GMS polar angle 9*, the resolution function 
is evaluated in separate bins of cos(0*). A simultaneous binned maximum likelihood 
£t is performed in each bin for all three channels, separately for SVDl and SVD2. 
The £t results are shown in Figure [T] 

Obtaining measurements of pcp and for all bins, the results are combined with 
least squares £t to constant to obtain the hnal result. 

The fit is tested on Monte Garlo simulation equivalent to six times the data statistics. 
A linearity test shows that no bias is present. 

2.2 Systematics and Final Result 

The estimated systematic uncertainties are listed in Table The main contribu¬ 
tion arises from the misalignment of the SVD and is estimated using Monte Garlo 
simulations for different misalignments. The uncertainty due to the position of the 
mass window is estimated by varying the position of the window. The uncertainty on 
background comprises a contribution due to statistical fluctuations and a contribution 
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Figure 1: Proper decay time distributions summed over cos{6*) bins and both SVD 
configurations with the sum of fitted functions superimposed. Shown as error bars 
are the distributions of events in the m{D^) signal region while the shaded area 
represents background contributions as obtained from sidebands. The corresponding 
pull is shown below each plot. 
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Source 

^Ycp (%) 

AA., (%) 

SVD misalignment 

0.060 

0.041 

Mass window position 

0.007 

0.009 

Background 

0.059 

0.050 

Resolution function 

0.030 

0.002 

Binning 

0.021 

0.066 

Total 

0.092 

0.066 


Table 1: Systematic uncertainties for the time-dependent analysis of —)■ 

K+K-.K-Ti+^T^+n-. 


that arises from modelling the distribution, which is estimated from MC simulation. 
The two contributions are added in quadrature. Systematics due to the resolution 
function are estimated using alternative parameterisations. Systematics due to bin¬ 
ning are estimated by varying the number of bins in cos(0*) and t. All individual 
systematic uncertainties are added in quadrature to obtain the overall uncertainty. 

The hnal result for the difference in the effective lifetime of mesons decaying 
to CP-even eigenstates iP+iP”, tt+tt” and the flavour eigenstate K~n~^ is 

ycp = [1.11 ± 0.22(stat.) ± 0.09(syst.)]% , 

which corresponds to a 4.7 a signihcance when statistical and systematic error are 
added in quadrature. Simultaneously, the CP asymmetry was measured, obtaining 

Ar = [-0.03 ± 0.20(stat.) ± 0.07(syst.)]% , (4) 

which is consistent with no CP violation. 


3 Time-integrated analysis of > K^K , tt+tt 

While the time-dependent analysis only measures indirect CP violation, it is possible 
to measure both direct and indirect CPV performing a time-integrated analysis of 
decay rates of neutral charm mesons to CP-even hnal states K~^K~ and 7r’''7r“. The 
presented results are an update of a 2008 analysis by the Belle Collaboration [5] 
using the full Belle data sample. The asymmetry of the time-integrated decay rates 
of 71° —)■ / and 71° —)■ / is dehned as 


P(7}0 ^ /) - r(71° ^ /) 

P(710^/) + P(:D°^7) 


(5) 
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where T is the partial decay width. The flavour of the charm meson is tagged via the 
charge of the slow pion from the decay D*~^ —)■ However, the experimentally 

measured quantity is 


N{do ^ /) - iV(D° ^ /) 
N{D^^f) + N(D°^7) 


( 6 ) 


where N is the number of events from a certain decay. This quantity, called raw 
asymmetry, comprises besides the physical CP asymmetry also the production and 
detector-induced asymmetry: Araw = Acp + ApB + A^±. The production asymmetry 
ApB is a forward-backward asymmetry in production of D*^ and D*~ and arises from 
7 — interference and higher order QED effects in e+e“ —)■ cc. It is assumed to be 
the same for all charm mesons. The term A^± is a detector-induced asymmetry that 
arises from different reconstruction efficiencies for positively and negatively charged 
particles. Since the decays in question are self-conjugated, the only charged 
particle that contributes to the asymmetry A^± is the slow pion, with A^± hence 
becoming A^{Tif). This term is determined using tagged and non-tagged decays 
—)■ K~7r~^, where it holds 


Auntag = ApB^D^) + Acp{Kn) + H,(i^7r) , (7) 

Atag = Afb{D*+) + Acp{K7i) + A,{K7l) + ) . (8) 


Since the forward-backward asymmetry is assumed to be the same for all charm 
mesons, it follows that = Atag — Auntag- Because A^{Kti) and A^iTif) are 

functions of the corresponding phase spaces in the laboratory frame, He(7r^) is cor¬ 
rected for separately in bins of and 6^,^. 

The forward-backward asymmetry is odd function of 9* and is corrected for using 

^CP = ilA“"(cosr) + /l“"(-cosr)| , (9) 

Afb = il9l“-(cos«-) - AZ:(- cos«-)| , (10) 

where A^^T, is the raw asymmetry after A^{Tif) correction. 


3.1 Event Selection 

The vertex fit and Pcms(D*’'') cuts are the same as stated in the previous Section. 
Other selection criteria are optimised so that the error on the asymmetry is minimal. 
The windows in m{D^) (q) are ±17.8 (1.00) MeV for K^K~, ±17.8 (1.85) MeV for 
K~tt^ and ±17.2 (0.90) MeV for 7r’''7r“. 

The signal yield is extracted via background subtraction. Background in the signal 
window is estimated based on two symmetrical sidebands that together amount to 
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Figure 2: CP-violating asymmetries in K^K and 7r'''7r final states, and forward- 
backward asymmetries for the same two final states. 


the width of the signal window, with their position starting at ±20 MeV from the 
nominal mass. Additionally, the range is limited to |cos(6'*)| < 0.8 in order to 
decrease Ae(7r^)-related systematics. 

The physical asymmetry Aqp is extracted in bins of cos(0*), with the hnal result 
obtained via a £t to constant. The corresponding plots are shown in Figure 

3.2 Systematics and Final Results 

The estimated systematic uncertainties are listed in Table The uncertainty cor¬ 
responding to the signal counting method comprises two contributions, summed in 
quadrature: possible small differences in the signal shapes of and D (uncertainty 
estimated from the tagged K~tt^ sample) and in background between the signal win¬ 
dow and sidebands. The uncertainty corresponding to sideband selection is estimated 
by varying the position of the sidebands. The uncertainty of the Ae(7r^) correction 
arises from statistics of the K~tt^ sample and binning (estimated from using different 
binnings and requiring different minimal statistics per bin). The extraction method 
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Ai-Acp (%) 

Signal counting method 

0.055 

0.023 

0.037 

ns correction 

0.065 

0.067 

0.014 

Aqp extraction method 

0.06 

0.050 

0.051 

Total 

0.085 

0.087 

0.064 


Table 2: Systematic uncertainties for the time-integrated analysis of —)■ 

K+K-,7l+7l-. 

of Acp contributes a systematic uncertainty due to binning and two different config¬ 
urations of SVD. 

The hnal results for the CP asymmetries are 

= [-0.32 ± 0.21(stat.) ± 0.09(syst.)]% , (11) 

A^Jp = [ 0.55 ± 0.36(stat.) ± 0.09(syst.)]% , (12) 

AAcp = [-0.87 ± 0.41(stat.) ± 0.06(syst.)]% , (13) 

and are consistent with no CPV. 

4 Time-integrated analysis of 

We report also a time-integrated analysis of CP violation in decays —)■ [B], 

performed on 966 fb“^ of Belle data. The analysis is largely similar to the time- 
integrated analysis described in the previous Section. 

4.1 Event Selection 

The is required to originate from the decay D*^ —)■ The same Pcms(D*’'') 

cut as before is applied. The optimisation of selection criteria is done in order to 
minimise the error on the raw asymmetry. The obtained range of m{D^) is (1.758, 
1.930) GeV and the q range is (0.14, 0.16) GeV. 

The corrections for production and detector-induced asymmetry is performed as pre¬ 
viously described. The signal yield is extracted via a simultaneous £t in Am of 
and D samples. The fit is performed in bins of (cos(6**),p^'®, cos{6'"^)), and the hnal 
result obtained via a £t to constant on values in bins of cos(6'*). The procedure is 
tested and confirmed on Monte Carlo simulations. The corresponding plots are shown 
in Figure 
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Figure 3: Projection of the simultaneous fit in Am to the (left) and (middle), 
and Acp in bins of cos(6'*) (right). 


4.2 Systematics and Final Result 

The estimated systematic uncertainties are listed in Table The hnal result for the 
CP asymmetry in decays —)• is 

Acp = [-0.03 ± 0.64(stat.) ± 0.10(syst.)]% , (14) 

which is consistent with no CP violation. It is the most accurate measurement to 
date. 


signal shape 

0.03 

Tis correction 

0.07 

Acp extraction method 

0.07 

Total 

0.10 


Table 3: Systematic uncertainties for the time-integrated analysis of D'^ —)■ 
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